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We present the results of a search for low energy Ve from the Sun using 1496 days of data from 
Super-Kamiokande-I. We observe no significant excess of events and set an upper limit for the 
conversion probability to Ve of the *B solar neutrino. This conversion limit is 0.8% (90% C.L.) of 
the standard solar model's neutrino flux for total energy = 8 MeV - 20 MeV. We also set a flux 
limit for monochromatic for Ea^ — lOMeV - 17MeV. 

PACS numbers: 14.60.Pq,26.65.+t,96.40.Tv,95.85.Ry 



Solar neutrino measurements at Super-Kamiokande T] problem is explained by the transformation of electron 
and SNO 2J have established that the solar neutrino neutrinos to other active neutrinos. The mechanism for 
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this transformation is generally assumed to be via neu- 
trino flavor oscillations from v^, to some superposition of 
v^^ and v-y. However, measurements reported so far do not 
rule out the possibility of spin flavor precession(SFP) in 
which some of the transform to antiparticles (i?^, Vr). 
In the so-called "hybrid models" [3] , SFP and oscillation 
can transform solar neutrinos to if the neutrino is Ma- 
jorana, it has a large magnetic moment, and the Sun has 
a large magnetic fleld. If the neutrino has a magnetic 
moment, there are two possibilities: (1) the neutrino is 
a Dirac particle ; (2) it is a Majorana particle. In the 
Dirac neutrino case, v]^ changes to by the spin mag- 
netic moment transition. The is a sterile neutrino. 
On the other hand, in the Majorana neutrino scenario, 
SFP causes Ve '^fi,T- Neutrino oscillation then yields 
j^fi.T — > f'e- Solar De could also originate from neutrino 
decay In this paper, we present a search for from 
the Sun. 

The inverse beta decay process, i>e+p — > n + e^ , is pre- 
dominant for De interactions in Super-Kamiokande (SK). 
The cross section for this process is two orders of magni- 
tude greater than that for elastic scattering, and there- 
fore SK has good sensitivity for the detection of solar Pg. 
The positron energy is related to the neutrino energy by 
Ef,+ « Ep^ — 1.3 MeV. The positron angular distribution 
relative to the incident direction is nearly flat with 
a small energy dependent slope |^, which is in contrast 
to the sharply forward peaked elastic scattering distribu- 
tion. The difference between these distributions can be 
used to separate solar neutrino events from events. 

Super-Kamiokande is a 22.5 kton fiducial volume wa- 
ter Cherenkov detector, located in the Kamioka mine in 
Gifu, Japan. The data used for the search were collected 
in 1496 live days between May 31, 1996 and July 15, 
2001. A detailed description of SK can be found else- 
where ■ Dominant backgrounds to the solar neutrino 
signal are ^^^Rn in the water, external gamma rays and 
muon- induced spallation products. Background reduc- 
tion is carried out in the following steps: first reduction, 
spallation cut, second reduction, and external 7-ray cut. 
The first reduction removes events from electronic noise 
and other non-physical sources, and events with poorly 
reconstructed vertices. The spallation cut removes events 
due to radio-isotopes (X) produced by cosmic ray muon 
interactions with water: /i -|-^^ O ^ fi + X . These radio- 
isotopes are called "spallation products." The spallation 
products emit beta and gamma rays and have lifetimes 
ranging from 0.001 to 14 sec. We cut these events using 
likelihood functions based on time, position, and muon 
pulse height. The time and position likelihood functions 
are measures of the proximity of a candidate event to 
a muon track, while the pulse height likelihood function 
measures the likelihood that a muon produced a shower. 
These three likelihood functions are used together to dis- 
criminate against spallation events 0. The second re- 
duction removes events with poor vertex fit quality and 
diffuse Cherenkov ring patterns, both characteristics of 
low-energy background events. The external 7-ray cut 



removes events due to 7-rays from the surrounding rock, 
photomultipliers(PMTs), etc.. Fig. ^ shows the energy 
spectrum after each reduction step. 
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FIG. 1: Energy spectrum after each reduction step. The solid 
curve shows the expected positron spectrum, after all cuts, 
assuming all *B solar neutrinos convert to Ue ■ The horizontal 
axis shows the reconstructed total energy. 

At SK, a positron from inverse beta decay is indistin- 
guishable from an electron or a gamma ray because the 
delayed 2.2 MeV gamma ray from n-l-p— >d-|-7is be- 
low the detector's energy threshold. In order to remove 
elastic scattering events due to solar neutrinos, we cut 
events with cos^sun > 0.5, where 9sun is the event direc- 
tion with respect to the direction from the Sun. The re- 
gion cos^sun < 0.5 would be occupied by solar events, 
in addition to events due to known background sources 
which could not be removed by the standard data reduc- 
tion. For E < 8 MeV, most background events are due to 
radioactivity in the detector materials (such as ^^^Rn). 
Spallation accounts for a small fraction of background 
events in this region. In contrast, for E > 8 MeV, most 
background events are produced by spallation. 

The spallation cut used in the data reduction efficiently 
removes short-lifetime spallation products. This cut also 
removes ^90% of long-lifetime products such as y^N (ri 

— 7.1 sec) and 4-^Be (ri = 13.8 sec). Event by event 
removal of the remaining ~ 10% of these events is im- 
practical because this introduces large dead time. How- 
ever, we can estimate the contribution of these events to 
the post-reduction data sample using a statistical sub- 
traction technique. First, we made a time distribution of 
muon events preceding each low energy event by up to 
200 seconds (Fig.El^A)). Since the average muon rate at 
SK is w 2.5 Hz, there are, on average, « 500 events for 
each low energy event. If the low energy event is due to 
a long lifetime spallation product, its event time will be 
correlated with one of the ~500 preceding muon events. 
If this is not the case, then its event time will be un- 
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correlated with all of the muon events. To estimate the 
number of /i responsible for spallation events, we have 
to subtract the number of /i which did not make spal- 
lation events from the total number of fi. In order to 
perform this subtraction, we made a sample of simulated 
events distributed randomly in space and time. We ap- 
plied the spallation cut to this sample as in the actual 
data sample in order to account for biases introduced 
by this cut. The muon time distribution of the random 
sample is shown in Fig. ll^B). The dip near Delta-T — 
is due to the accidental loss of events by the spalla- 
tion cut. To estimate the number of muons which made 
spallation products, distribution (B) with suitable nor- 
malization is subtracted from distribution (A); the result 
of this subtraction is shown in Fig.l^C). The number of 
muon events in the delta-T =100 sec - 200 sec region 
is used to calculate the normalization factor because the 
contamination from muons which make spallation prod- 
ucts is negligible in this region. The number of spallation 
events is obtained as 
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Nq'^^qs^^ is the number of muon events within 50 seconds 
preceding the observed events, while -/Vq^^q"™ is the cor- 
responding number for random events. N°QQ'^200sec 
■^iSo-200sec ^-^s similarly defined, but with a timing win- 
dow of 100 to 200 seconds preceding the events. For 
8.0-20.0 MeV, and cos^sun < 0.5, the number of spalla- 
tion background events obtained by this method is (2.77 
± 0.20) X 10"*. The number of observed candidate 
events is 29781, so the ratio of spallation events to ob- 
served events is (93±7)%. The spallation contamination 
in each energy bin is shown in Fig. |31 

The energy spectrum of the solar is not known be- 
cause the mechanism for D,, creation is not known. Even 
if one assumes the SFP-oscillation hybrid model, the en- 
ergy spectrum depends on x Bso/ar, Am^ and sin^(26'), 
none of which are known precisely, if at all. In order to 
deal with this ambiguity, we have chosen two spectrum 
models: the ®B neutrino spectrum 7] and monochro- 
matic spectrum (spectrum independent analysis). 

For the ^B spectrum dependent analysis, we obtain 
an upper limit on the solar D,, flux by comparing the 
observed number of events outside of the elastic scatter- 
ing peak (cos^sun < 0.5) with the expected number of 
Ve events assuming that all ^B neutrinos convert to Vf.. 
The expected number is obtained by Monte Carlo sim- 
ulation of solar Df. interaction with the detector. The 
cos ^sun dependence was simulated , and the effect of this 
dependence on the Vf, efficiency is taken into account. 
The standard solar model(SSM) ®B neutrino flux was as- 
sumed (5.05 xlO^ /cm^/sec) |g] . Through the rest of 
this paper, electron neutrino spectrum and flux refer to 
the unoscillated quantities at the Sun. The solid lines in 
Fig. 4 show 90% C.L. limits on the Vf. flux before statis- 
tical spallation subtraction. The dashed lines show the 
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FIG. 2: (A): ^ delta-T distribution before observed events 
(B): Before random events (C); The delta-T distribution of 
events caused by spallation products obtained as (A) — scale 
factor x (B). 
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FIG. 3: Spallation contamination in each energy bin. The 
horizontal axis shows the total energy and the vertical axis 
shows the ratio of spallation events to observed events. 



limits after statistical subtraction (only for E > 8 MeV) . 
By combining the statistics for 8 MeV < E < 20 MeV, 
we obtained a global upper limit of 0.8% of the SSM 
neutrino flux. 

Some authors have indicated that the positron angu- 
lar distribution may be useful for the search for in 
the SK data (e.g. 0, 0|). cos^sun is distributed as 
f(cos6'sun) = 0.5 X (1 -I- a X cos6'sun), where a is a mono- 
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tonically increasing function of neutrino energy (except 
near threshold), and a < for < 13 MeV and > 
above this 0. The angular information is useful for the 
Df, search at the lowest neutrino energies where f(cos 0sun) 
has sufficient slope and the event statistics are large. Ve 
events with the predicted cos 0sun distribution were input 
to a detector simulator to obtain the expected positron 
angular distribution. The resulting distribution has the 
same form as above. The fit value of a is —0.076 at i? = 
5 - 6 MeV, 0.107 at £; = 12 - 20 MeV, and crosses at 
- 9 MeV. 

Solar neutrino elastic scattering is one of the back- 
grounds for this analysis. Almost all such events have 
cos^sun > 0.5, so events with cos^sun > 0.5 are cut. We 
also subtract the small amount of spill-over into cos ^sun 
< 0.5 using Monte Carlo simulation (--5% for 5-20MeV). 
Another background is due to ^^0{Vf.;eY^¥ 11]. There 
is only a small number of events from this source (0.03% 

2%, depending on energy), but electrons from this pro- 
cess, like the low-energy j>e, have negative slope in their 
angular distribution. So they are subtracted from the 
data. The flux is taken as the charged current flux 
value from SNO, 1.76 x 10^ /cmVsec |l^ . 

A Ve flux upper limit is obtained using a probability 
test with the slope of the cos 0sun distribution serving 
as a constraint. This test is based on a test with 
X^((5, /3,7) defined for each energy as follows: 



the spallation background subtratction is not applied in 
this analysis for two reasons. First, for E < 8 MeV, spal- 
lation subtraction is ineffective because spallation events 
form a small subset of the total background. Second, for 
E > 8 MeV, there are insufhcient statistics after spal- 
lation subtraction to perform an angular analysis of the 
data. 
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i is the index for the cos6'sun bins (cos^sun < 0.5, Ncos= 
30), Xi is (cos0sun)z, Nf"-*'^ is the number of observed 
data events, a?*"*- is the statistical error of the observed 
data, Nf is the expected number of elastic scattering 
events, ^ is the expected number of events from the 
^*0(z^e;e)^*F reaction, N^^ is the number of events 
assuming all SSM convert to Ve (the number in each 
bin i depends on the slope a), and nf^ is the shape for 
all other background events that are almost uncorrelated 
in direction with the Sun (this background is essentially 
flat). Nf and N-^^ are both < 2% of iV^'^*^ and the 
systematic errors of these terms are negligible. cr^^** (= 
0.5%) is the systematic error of the slope of the back- 
ground shape and 7 is the parameter that takes this into 
account. /3 parameterizes the amount of such background 
events. We divided the parameter space for S into a grid, 
and minimized with respect to /3 and 7 at each grid 
point. The resulting 7 and Xmin indicated good fits to 
the data, as a function S obtained in this way is input 
to a probability function. From this analysis, we set a 
90% C.L. upper limit for each energy bin. The dotted 
lines in Fig. ^ show the result. It should be noted that 



FIG. 4: Summary of limits. The horizontal axis shows 
total positron energy and the vertical axis shows the 90% 
C.L. z/e rate normalized to the SSM Ve rate. The solid lines 
show the 90% C.L. limit ratio. The dashed lines show the 
hmit after statistical subtraction of the spallation background. 
The dotted lines show the result from the angular distribution 
analysis. 

The analysis described so far assume that the orig- 
inate from solar neutrinos. We also generalized our 
search by assuming a monochromatic source at var- 
ious energies and set conservative flux upper limits. 
The interaction of such with the detector was sim- 
ulated, and standard data reduction cuts were applied. 
The positron spectrum is well described by a Gaussian. 
We then counted the number of events in the data in the 
± 1(7 range of this Gaussian. We took this number to 
be the number of events due to monochromatic De and 
obtained an upper limit. This upper limit is very con- 
servative because we do not take account of the large 
spill-over from lower energy bins that is implied by the 
sharply falling spectrum seen in the data. We also ob- 
tained limits after statistical subtraction of long lifetime 
spallation events. The 90% C.L. limits are shown in Fi^ 

In summary, a search for flux from the Sun was per- 
formed using all 1496 live days of solar neutrino data from 
Super-Kamiokande-I. Using the ®B and monochromatic 
energy spectra, 90% C.L. upper limits were set for the De 
flux. For the ®B spectrum dependent analysis, the upper 
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FIG. 5: Pe flux 90% C.L. upper limit for each monochromatic 
v^. The horizontal axis shows neutrino energy and the vertical 
axis shows the flux limit. The circles show the limits before 
spallation subtraction while the stars show the limits after 
subtraction. The two highest-energy bins have insufficient 
number of events for statistical subtraction. 



limit to the flux was 0.8% of the SSM Vf. flux prediction 
for E = 8.0-20.0MeV. This can be compared with the 
Kamiokande result of 4.5% For v^, fluxes with var- 
ious monochromatic energies, the resulting upper limits 
are shown in Fig. |5| 
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